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^ ; Abstract 

We discuss modeling of generalized parton distributions (GPDs), their access from 
present experiments, and the phenomenological potential of an electron-ion col- 
' lider. In particular, we present a comparison of phenomenological models of GPD 

Q>^ ■ H, extracted from hard exclusive meson and photon production. Specific em- 

OO , phasis is given to the utilization of evolution effects at moderate xbj in a future 

' high-luminosity experiment within a larger lever arm. 

vn ■ 
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1 Introduction 

■ 

^ ■ Generalized parton distributions (GPDs), introduced some time ago [HO [3], have received 

, much attention from both the theoretical and the experimental side. This was triggered 

by the hope to solve the 'spin puzzle', referring to the mismatch of quark contribution 
to proton spin, as extracted from polarized deep inelastic scattering, and as given by the 
constituent quark model. We view the 'spin puzzle' first and foremost as a quest to quantify 
the partonic structure of the nucleon in terms of quark and gluon angular momenta, where 
an appropriate decomposition of the nucleon spin in terms of energy-momentum tensor form 
factors A'^{Q'^) and B'^{Q'^) has been suggested by X. Ji [4j: 

,=uA- (1) 
J\Q') = 2 [^'(2') + B'(Q')] = y J [^'(^' V,t = 0, Q2) + E\x, r?, t = 0, . 

The quark and gluon contributions are given by the first moments of parity-even and target 
helicity (non-)conserved GPD H (E). 

Furthermore, it has been realized that GPDs allow for a three-dimensional imaging 
of nucleons and nuclei [5], providing, in the zero-skewness case (r/ = 0), a probabilistic 
interpretation in terms of partonic degrees of freedom [6j. By definition GPDs are linked to 
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parton distribution functions (PDFs) and elastic form factors. In phenomenology they are 
used for modeling elastic form factors and the description of hard exclusive leptoproduction 
or even photoproduction. For these hard exclusive processes factorization theorems have 
been proven in the collinear framework at twist-two level [71[8]. In fact, GPDs build up a 
whole framework for description of hadron structure [HI [TU], with the 'spin puzzle' being 
just one interesting aspect. 

Much effort to measure hard exclusive processes has been spent in the last decade by 
the HI and ZEUS collaborations (DESY) in the small xbj region and at the fixed target 
experiments HERMES (DESY), CLAS (JLAB), and Hah A (JLAB) in the moderate xbj 
region. Thereby, deeply virtual Compton scattering (DVCS) off nucleon is considered as 
the theoretically cleanest process offering access to GPDs. Its amplitude can be parame- 
terized by twelve Compton form factors (CFFs) jllj , which are given in terms of twist-two 
(including gluon transversity) and -three GPDs. E.g., at leading order (LO) parity-even 
twist-two CFFs, Ti and £, can be expressed through quark GPDs H and E and they take 
the form: 

{ ^ } (xBj ,t,Q') f dx ^T^lr— { f } ^ = Q') . (2) 

Here both quark and anti-quark GPDs might be defined in the region x € [— ^, 1]. Sim- 
ilar expressions can be written for twist-two parity-odd CFFs % and £, while for other 
CFFs they are a bit more intricate [H]. The Bjorken variable xbj might be set equal to 
2^/(1 + ^). Analogous formulae hold for the LO description of 7*A^ — )• MN transition 
form factors (TFFs), measurable in deeply virtual electroproduction of mesons (DVEM). 
Here, in addition to GPDs, the non-pertubative meson distribution amplitude enters, which 
describes the transition of a quark-antiquark state into the final meson. This induces an 
additional uncertainty in the GPD phenomenology. 

Let us shortly clarify which GPD information can be extracted from experimental mea- 
surements. Neglecting radiative and higher twist-contributions, one might view the GPD 
on the rj = X cross-over line as a "spectral function", which provides also the real part of 
the CFF via a "dispersion relation" [HI [H ESI [I7j : 

QmT{xBi,t,Q^) ^= 7rF{ttt,Q^), F = {H,E,H,E}, (3) 
3^e|^}(xBj,t,Q2) = F\l^'dx^^!^^^yx,x,t,Q')±V{t,Q^). (4) 

The GPD support properties ensure that Eqs. ([3]) and @ are in one-to-one correspondence 
to the perturbative formula ([2]), where the subtraction constant V, related in a specific 
GPD representation to the so-called D-teim [12], can be calculated from either H or E. 
However, we note that the "dispersion relation" (j4]) differs from the physical one by the 
support property of the spectral functioij^]- To pin down the GPD in the outer region 
y > 7] = X, one might employ evolution, e.g., in the non-singlet case the change of the GPD 

^The physical or hadronic dispersion relation possesses a t/Q'^ dependent threshold that approaches one 
in the (generalized) Bjorken limit —>■ oo, see e.g., Ref. |14) . It is obvious that the twist expansion of 
the DVCS amplitude induces this threshold artifact on partonic level, related to the intricate problem of 
higher twist contributions [TB], where even the choice of partonic momentum fraction or scaling variable 
is nontrivial. We emphasize that for massless pions and within the setting XBj = 2^/(1 + the partonic 
"dispersion relation" Q yields the hadronic one, written in terms of XBj within the support < xbj < 1, 
where, however, both the spectral functions and integral kernels differs by i/Q^-dependent terms. 
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on the cross-over line is governed by (the equation in the whole outer region is needed) 

/i^-^F(x,x,t,/u^) = / —V{l,y/x,asiiM))F{y,x,t,fi'^). (5) 
a/i X 

Here, the kernel might be written to LO accuracy as [l]: 

Vil,z>l,a,) = ^!^j^ + J'dz' ^6(1 -z) + h{l - z) I + Oial) . (6) 

Unfortunately, a large enough range is not available in fixed target experiments. Hence, 
we must conclude that in such measurements essentially only the GPD on the cross-over 
line (thanks to also outside of the experimentally accessible part of this line flTj) and 
the subtraction constant V can be accessed. Moments, such as those entering the spin sum 
rule ([1]), can only be obtained from a GPD model, fitted to data, or more generally with 
help of some 'holographic' mapping |17j : 

{F{x,r] = 0,1,0"), F{x,r] = x,t,Q^)} =^ F{x,7],t, Q^) . (7) 

Here, F'^{x,rj = 0, Q^) are constrained from form factor measurements and, additionally, 
GPDs [H^) by (un)polarized phenomenological PDFs. Of course, a given 'holographic' 
mapping holds only for a specific class of GPD models. 



2 GPD modeling 

The implementation of radiative corrections, even including LO evolution ([5]), requires to 
model CFFs or TFFs in terms of GPDs. This can be done in different representations, 
which should be finally considered as equivalent. However, for a specific purpose a particular 
representation may be more suitable than the others. 

First, GPDs might be defined as Radon transform of double distributions (DDs) [H 118]: 

F{x,ri,t,i?)= 1 dy I dz {1 - xf5{x - y - zr])f{y, z,t, fi") , (8) 
Jo J-i+y 

where integer p G {0, 1}. In this representation polynomiality, however, not positivity con- 
straints are explicitly implemented. Moreover, with the right choic^ for p, the polynomiality 
of x-moments can be completed to the required order in rj. In the central, — r/ < x < r], 
and outer, r] < x < 1, region the GPD can be interpreted as the probability amplitude of a 
t-channel meson-like and s-channel parton exchange, respectively. Mathematically, F is a 
twofold image of the DD /, where the central and the outer region can be mapped to each 
other [261 [T^ [T7| . The potential ambiguity, given by a term that lives only in the central 
region, is removed by requiring analyticity [IH [T7] . 

Popular GPD models are based on Radyushkin's DD ansatz (RDDA) [T8] for t = 0, 
where the DD factorizes into the PDF analogue f{y) and a normalized profile function 
n(z). The GPD on the cross-over line is then given as 

/'^/(£ll^)n(.), (9) 

J_i I — xz \ I — xz J 

^Note that for the GPD E the factor (1 — x) is suggested by a spectator quark model analysis |19) : 
however, it might be replaced by a more general first order polynomial. For instance, for a spin-zero target 
the choice x looks rather natural [24], which has been recently discussed in detail in Ref. |25) . 
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which is a hnear integral equation of the first kind within the kernel — xz). 

Knowing the GPD at rj = 0, i.e., /(y), and on the cross-over line, allows to determine the 
profile function and so to reconstruct the entire gpeI, giving example of the 'holographic' 
mapping ([7]). 

On the first glance a GPD in the outer region can be straightforwardly represented by 
an overlap of light-cone wave functions (LCWF) [271 E5], which guarantees that positivity 
constraints are implemented. In simple models one might even reduce the number of non- 
perturbative functions, e.g., in a spectator diquark model, one only deals with one effective 
scalar LCWF for each struck quark species. This predicts for each parton species four 
chiral even and four chiral odd GPDs. Also one might utilize the overlap representation 
to evaluate both GPDs and transverse momentum dependent parton distributions (TMDs) 
from a given LCWF model. Unfortunately, there is a drawback. In the central region the 
GPD possesses an overlap representation in which the parton number is not conserved, 
and where the LCWFs are dynamically tied to those used in the outer region. A closer 
look reveals that Lorentz covariance already ties the momentum fraction and transverse 
momentum dependence of a LCWF [19j, see also the work in Refs. [201 [211 [22l [23] . Hence, a 
overlap representation is only usable if the LCWFs respect Lorentz symmetry, which would 
allow to restore the GPD in the central region [19] ■ 

Strictly spoken, positivity constraints for GPDs are only valid at LO, since they can be 
violated by the factorization scheme ambiguity. Nevertheless, it would be desired to impose 
them on GPD models. One might follow the suggestion of [29j and model GPDs as an 
integral transform of (triangle) Feynman diagrams, i.e., spectator quark models. A specific 
integral transformation, namely, a convolution with a spectator mass spectral function, can 
be used to include Regge behavior from the s-channel view [30j. Such dynamical models 
provide also effective LCWFs or TMDs; however, simplicity is lost. In particular, PDF 
and form factor constraints cannot be implemented, i.e., one has to pin down such models 
within global fitting. 

At present we neglect positivity constraints and we model GPDs in the most convenient 
manner by means of a conformal SL(2,M) partial wave expansion, which might be written 
as a Mellin-Barnes integral |26j 

2 y^-ioo sm(7rj) 

Here, pj{x,r]) are the partial waves, given in terms of associated Legendre functions of 
the first and second kind, and the integral conformal GPD moments Fj{r],t, fj?) are even 
polynomials in ry of order j or j -|- 1. We note that various other representations are based 
on the SL(2,M) partial wave expansion, see e.g., Ref. [2B]. In particular, the so-called "dual" 
parametrization [31], initiated in Ref. [32], has been confronted with the RDDA [331 [M]. 

In the Mellin-Barnes representation the CFFs possess a rather convenient form, e.g., 
Eq. (|2]) is rewritten as [HSlfH] 

''An example is provided by f{x) oc x~°'{l — x)^, which yields after some redefinitions the integral kernel 
fe = (1 — xz)°'~'^~^. The solution is then obtained in Mellin space and can be given as a convolution integral. 
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Figure 1: A simple GPD model (long dashed), based on the ansatz ([T2|) . versus Alekhins LO PDF 
parameterization [39] (grayed area) [left panel] and Kelly's [40] (dotted) [Sachs (short dashed)] form 
factor parameterization [right panel]. 



This integral is numerically implemented in an efficient routine in two different factorization 
schemes, including the standard minimal subtraction (MS) one at next-to-leading order 
(NLO) accuracy. Further advantages of this representation are: 

• The conformal moments evolve autonomously at LO. 

• One can employ conformal symmetry to obtain next-to-next-to-leading order (NNLO) 
corrections to the DVCS amplitude [35] [36], 

• PDF and form factor constraints can be straightforwardly implemented. Namely, 
Fj{7] = 0, t = 0, fi^) are the Mellin moments of PDFs, Fj=Q are partonic contributions 
to elastic form factors, Hj=i and Ej=i are the energy-momentum tensor form factors, 
and for general j one immediately makes contact to lattice measurements. 

Let us now illuminate the GPD model aspect based on generic arguments and a simple 
ansatz, e.g., for valence-like contributions at an input scale fiQ-. 



r(l - a{t) + j)r(l +p-a + j)T{2 - a + /3) 



r(i- 
(1 



a)T{l+p- 

r(2 

-h) + h^ 



a{t)+j)T(2-a 
-a + p + 5p)T{2- 



-j + P) 
a + j + P) 



r(2 - a + i + /3 + 5/3)r(2 - q + /3) 



(12) 



Here, n = 2(1) for u (d) quarks, a{t) = a + a't is the leading Regge trajectory, p determines 
the large —t behavior, /3 and 6^ the large j behavior, and /i is a phenomenological parameter. 
The PDFs are then given by an inverse Mellin transform and read: 



n 



r(2-a + /3) 



r(i- 
(1 



"°(i - 

a)r(l + /3) ^ 



x)^ 

-a + ^ + SP) 



(13) 



r(2 - a + /3)r(l + /? + 5/3) 



(1-x) 



5/3 



In the case that large x counting rules [37| would not be spoiled by non-leading terms in a 
1 — X expansion, h x n might be interpreted as the amount of a quark which has opposite 
helicity to the longitudinally polarized proton. That such a model provides reasonable 
results has been argued for the iso-triplet part of H and E [38]. In Fig.[T]we illustrate this for 
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the valence part of the proton GPD H = (4/9)/?" + {1/9)H and the electromagnetic 
form factor Ff, where the model result is shown as long dashed curve. To adopt to Alekhins 
LO PDF parameterization [35] we choose the Regge intercept a = 0.43, /3 = 3.2, 6f3 = 2.2, 
h = —1. For the form factor we take the slope parameter a/ = 0.85 and with the choice 
p = 2.12 the outcome is hardly distinguishable from Kelly's parameterization (dotted) |40) . 
Note that /3, 5f3, and p only differ slightly from the canonical values 3, 2, and 2, respectively, 
and that a{t) = 0.43 + 0.85t is essentially the p/uj trajectory. Moreover, an inverse Mellin 
transform provides the t-dependent zero-skewness GPD, which has been also alternatively 
modeled in momentum fraction space |41l H2] . 

To parameterize the degrees of freedom that can be accessed in hard exclusive reactions, 
one might expand the conformal moments in terms of t-channel SO (3) partial waves [43] 
dj{r]), expressed by Wigner rotation matrices and normalized to dj{r] = 0) = 1. An effective 
GPD model at given input scale Qq is provided by taking into account three partial waves, 

Fj{v,t) = dMft'it) + rfd,.2^)ff'{t) + v%~4{r^)ff'{t) , (14) 

valid for integral j > 4. In the simplest version of such a model, one might introduce just 
two additional parameters by setting the non-leading partial wave amplitudes to: 

fr''ir],t)=Skf^+'(rj,t), k = 2,4,.--- (15) 

Such a model allows us to control the size of the GPD on the cross-over line and its Q^- 
evolution, see Figs. [3] and [TT| for small or moderate x values, respectively. A flexible 
parameterization of the skewness effect in the large x region requires to decorate the skew- 
ness parameters Sk with some j dependence and for more convenience one might replace 
Wigner's rotation matrices by some effective S0(3) partial waves. 



3 GPDs from hard exclusive measurements 

Based on the experimental data set from the collider experiments HI and ZEUS at DESY, 
the fixed target experiment HERMES at DESY, and the Hall A, CLAS, and Hah C ex- 
periments at JLAB, GPDs have been accessed from hard exclusive meson and photon elec- 
troproduction in the last few years. Favorably, DVCS enters as a subprocess into the hard 
photon electroproduction where its interference with the Bethe-Heitler bremsstrahlung pro- 
cess provides variety of handles on the real and imaginary part of twist-two and twist-three 
related CFFs |441lllj. However, switching from a proton to a neutron target allows only for 
a partial flavor separation, which is much more intricate than in deep inelastic scattering 
(DIS). On the other hand DVEM can be used as a flavor filter, however, here one expects 
that both radiative [451 H6l BTj and (non-factorizable) higher-twist contributions might be 
rather important. The onset of the collinear description remains here an issue, which should 
be phenomenologically explored. 

For the DVCS process, the collinear factorization approach has been employed in a spe- 
cific scheme up to NNLO in the small xbj region [351 |36l [Hj. It turns out that NLO cor- 
rections are moderate, while NNLO ones are becoming much smaller |14| . Experimentally, 
the unpolarized DVCS cross section has been provided by the HI and ZEUS collaborations 
\4:8\ H9l [50l [51] . In these collider kinematics the cross section is primarily given in terms of 
two CFFs, nandS: 



-{W,t,Q 



(^Bj,t,Q') 32 • (16) 
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Although at a fixed scale and to LO accuracy the CFFs are given by (dominant sea) quark 
GPDs, evolution will induce a gluonic contribution, too. Indeed, the experimental lever 
arm 3GeV^ ^ ^ 80GeV^ is sufficiently large to access the gluonic GPD. In a fitting 
procedure the Mellin-Barnes integral was utilized within a S0(3) partial wave ansatz for 
the conformal moments and good fits (x^/d.o.f. ~ 1) could be obtained at LO to NNLO 
accuracy, exemplifying that flexible GPD models were at hand. From such fits, one can 
then obtain the image of quark and gluon distributions. It is illustrated in Fig. [2] that in 
impact space, the (normalized) transverse profiles. 



p(b,x,Q^) = ^^ ^ ^ (17) 



determined for dipole and exponential t-dependence of H, mainly differ for distances larger 
than the disc radius of the proton, i.e., for b > 0.6 fm. Hence, the larger values of the 



transverse widths \/ (b"^) ~ 0.9 fm and \/ (6^) ~ 0.8 fm for the dipole ansatz arise from 

V sea V G 

the long-range tail of the profile function, see the solid curve. For an exponential ansatz we 



find slightly smaller values y {b"^) ~ 0.7 fm and y (6^) ~ 0.6 fm, where the gluonic one 

V SG& V ^^^^ Cj 

is compatible with the analysis of J/^ production [52]. Note that the model uncertainty 
in the extrapolation of the GPD to t = corresponds to the uncertainty in the long-range 
tail. Moreover, the model uncertainty of the extrapolation into the region —t>l GeV^ is 
essentially canceled in the profile (jl7p . normalized at 6 = 0. 

Also the model dependent extrapolation 77 — >■ has been employed above, which is 
controlled by the skewness effect. This effect might be quantified by the ratio 

^^ Hix,xt = OQ^) ^ g(x,Q2) = F(x,0,t = 0,Q2). (18) 
q{x, Q^) 

In a minimalist SO (3) partial wave model this ratio is at small x determined by the effective 
Regge intercept and one finds for sea quarks (gluons) the conformal, so-called Shuvaev, 
ratio ~ 1.6 (~ 1.1) [5l], see dotted curves in Fig. [3j These ratios arises in a large class of 
GPD models, including the conformal one, and were in the past widely misunderstood as a 
prediction that ties GPDs and PDFs in the small 2;-region. Contrarily, we have found to LO 
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Figure 3: The r ratio (|18l) of sea quark (left) and gluon (right) GPDs for a; = 10 and t = versus 
Q2 for Giro? model 58 (dashed) and two of our flexible GPD models: KM 10a (sohd) and KMlOb 
(dash-dotted). The conformal ratio is shown as dotted curve. 



accuracy that the r-ratio is for quarks approximately one and rather stable under evolution, 
while for gluons it is much smaller than one. This is illuminated in the left and right panels of 
Fig. El respectively, where we utilized two models KMlOa (solid) and KMlOb (dash-dotted) 
that take into account three S0(3) partial waves, see below. Thereby, the gluonic value 
r*^ < 1 ensures the stability of the sea quark r^-ratio under LO evolution, where, however, 
its precise value cannot be pinned down. Even a negative value at < 5 GeV^ (solid) 
is compatible with present DVCS data, which might be considered as a model artifact. 
We also note that at LO the gluonic GPD as the gluonic PDF are rather steep and that 
radiative corrections might provide a large GPD/PDF reparameterization effect, which will 
be in future studied in more detail. Our first successful LO description of DVCS within 
a flexible GPD model is in agreement with aligned-jet model considerations [55]. We also 
mention that an attempt has been undertaken to access the £ OFF from the beam charge 
asymmetry measurement [56], proportional to the combination Ke [Fi{t)T-L — j^F2{t)£~\ . 
Unfortunately, the size of experimental uncertainties does not allow to separate the Ti and 
£ contributions. 

An approach analogous to the one employed for DVCS in Ref. [H] is also suitable for 
LO and NLO analysis of DVEM. Hence, one can simultaneously make use of DVCS and 
DVEM measurements in a global fitting procedure, which is in progress. The hard exclusive 
vector meson production has been also extensively studied in both the color dipole and the 
hand-bag approach. In the latter approach [53 [58] a modified factorization scheme, used at 
LO, has been assumed in which the transverse size of the meson is taken into account. Note 
that in such an approach radiative corrections, mainly related to the produced meson state, 
are partially taken into account, while the GPD is treated in the collinear approach. The 
GPD model GK07 from Ref. [58j is based on RDDA, where the PDF parameterizations, 
including the evolution, are taken from the CTEQ6M NLO fit [59]. According to the 
authors, this skewing of Q^-evolved PDFs at NLO is consistent with the evolution of GPDs 
to LO accuracy. Consequently, at t = the H GPDs are already fixed and inherit the PDF 
properties. It is illustrated in Fig. H] (left) that for > 5 GeV the data are well described 
in this approach. In the large W region the cross section is gluon dominated; however, 
sea quarks play an important role, too (about 40% for production longitudinal cross 
section in HI and ZEUS kinematics). The dashed curve shows a GPD LO model prediction 
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Figure 4: Longitudinally total cross section of 7lP — >■ p^p (left) at = 4GeV^ and (jyp 
at Q2 = 3.8 GeV^ (right) versus W from HI [BOl EH [62] (filled squares), ZEUS [63l [64] (open 
squares) and [65] (open triangles), CLAS [66] (open circles), HERMES [HT] [68] (filled circles), and 
CORNELL [69] (filled triangles). Curves display predictions from the hand-bag approach with the 
GK07 model [5S] (solid) and the coUinear LO approach within the KMlOb model from a DVCS fit 
[70] (dash-dotted). (Figures are taken form Refs. [71j and l72i.) 



in the collinear factorization approach, where the gluons have been extracted from small 
XBj DVCS data via scaling evolution effects. In a LO description of DIS the gluon PDF 
is rather steep and although our skewness ratio is ~ 0.5, we certainly fail to describe the 
data. We emphasize that the skewness ratio for sea quarks in the GK07 parametrization 
lies below the conformal one and is even rather stable under evolution, cf. LO evolution 
examples within RDDA in Fig. 6 of Ref. |73] . These properties of our GPD models and 
partially of the GK07 in the quark sector are qualitatively different from the properties 
of the Durham GPD parameterization [74] . which, e.g., is employed in the diffractive J/^ 
electroproduction [75], and provides the conformal ratio, see also dotted curves in Fig. [3j 
In the resonance region the GK07 model does not describe the production data, but it 
does describe the (j) production. Whether the approach is not applicable in this region or 
the valence-like GPDs looks rather different from those obtained with the RDDA, is for us 
an open question. 

GPD studies were also performed for the DVCS process in the fixed target kinematics 
to LO accuracy. In this region, relying on the scaling hypothesis, one might directly ask for 
the value of the GPDs on their cross-over line. For instance, for valence quarks we use the 
generically motivated ansatz 

^...(....0 = 1^ f^V°"' fi^V f 1 - i^^V . (19) 

Here, the skewness ratio r = limx^oH{x,x)/H{x,0), a{t) = 0.43 -|- 0.85t/GeV^, b controls 
the X — > 1 limit, and the residual t-dependence, which we set to M™^ = 0.8 GeV, 

where q{x) = H(x,0) is a reference PDF, e.g., the LO parameterization of Alekhin [39j . 
The generic (— 1)~^ fall-off at large —t for generalized form factors is indirectly encoded 
in the Regge-trajectory and the residual t dependence, chosen by a monopole with an x- 
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dependent cut-off mass. The subtraction constant is normalized by d and controls the 
i-dependence: 



V{t) = d{l 



Mi) 



-2 



(20) 



In a first global fit ^3] to hard exclusive photon electroproduction off unpolarized pro- 
ton we took sea quark and gluon GPD models with two S0(3) partial waves at small x, 
reparameterized the outcome from HI and ZEUS DVCS fits at = 2 GeV^, and employed 
it in fits of fixed target data within the scaling hypothesis. To relate the CFFs with the 
observables we employed the BKM formulas within the 'hot-fix' convention |76j and 
used the Sachs parameterization for the electromagnetic form factors. Thereby, we utilized 
the "dispersion relation" (I3|4p . where the ansatz pOh specifies a valence- like GPD on the 
cross-over line. Besides the subtraction constant ()20p . we also included the parameter- free 
pion-pole model for the E GPD \77\ and parameterized the H GPD rather analogously to 
Eq. ([19]) with b = 3/2. For the fixed target fits we chose two data sets, resulting in two fits 
{KM09a and KM09h). The first set contains twist-two dominated (preliminary) beam spin 
asymmetry. 



^^^^'^^ {dxBjdtdQ^d4 
= 4)sin(,^) + .. 

and beam charge asymmetry, 




dxBjdtdQ 



(21) 



AbM = 



da^ 
^dxBjdtdQ'^ 



= 4°)+4)cos(<^) + ... 




(22) 



coefficients a'^I and A^j^l, respectively, from HERMES [751 [TS] and 12 beam spin asymmetry 

coefficients A^g, which we obtained by Fourier transform of selected CLAS |80j data with 
small —t. The second data set includes also Hall A measurements [81] for four different t 
values. In light of the discussion [82J of Hall A data, we projected on the first harmonic of 
a normalized beam spin sum 



^BS 



27r 



dw cos 




(23) 



where dw oc Vi{(p)V2{(p)d(f) includes the Bethe-Heitler propagators. We haven't used Hall 
A helicity-dependent cross sections (beam spin differences). The Hall A data, given at 
relatively large xbj = 0.36, can only be described in our model within an unexpectedly 
large value of H, which is not visible in single longitudinally target spin asymmetries at 
smaller values of XBj ■ Otherwise, our findings 
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Figure 5: Experimental measurements for fixed target kinematics (circles) labeled by data point 
number n: A'^^ (1-18), (19-36), (37-54) from Ref. (SB]; A^^^ (55-66) and E^^^'"" (67-70) 

are derived from Refs. [80 and [81 . Model results are pinned down by "dispersion- relation" fits 
KM09a without (squares, slightly shifted to the l.h.s.) and KM09b with (circles, slightly shifted 
to the r.h.s.) Hall A data [53], hybrid model fit KMlOb (triangles- up), and a hand-bag prediction 
GK07 irom hard vector meson production (triangles-down, slightly shifted to the r.h.s.) [55] . 



are compatible with our generic expectations: the skewness effect at small x should be small, 
i.e., r ~ 1, the subtraction constant should be negative [Ml El]) and, according to counting 
rules [85], b should be smaller than the corresponding f3 value of a PDF, see Refs. [53l [70]. 

To improve just described models, we now use a hybrid technique where sea quark and 
gluon GPDs are represented in terms of conformal moments, while, for convenience, the 
valence quarks are still modeled in momentum fraction space and within the "dispersion 
integral" approach. Also, the residue of the pion-pole contribution is now considered as a 
parameter, and the Hall A data forces a roughly three times larger value than expected 
from the model [77]. Optionally, we might also use the improved formulae from Ref. [87] . 
applicable for a longitudinally polarized target. The parameters, 

KMlOa: r™' = 0.88 , M™^ = 1.5 GeV , 6™' = 0.40 , d = -1.72 , Md = 2.0 GeV , 
KMlOh: r™' = 0.81 , M^^ = 0.8 GeV , 6™' = 0.77 , d = -5.43 , = 1.33 GeV , 

for the valence part of H GPD are qualitatively compatible with those from the pure KM09 
"dispersion relation" fits. 

We also did one fit where we directly used harmonics of beam spin sums and differences 
measured by Hall A (e.g. numerator of r.h.s. of ()23p ). obtaining parameters 

KMIO: r™^ = 0.62 , M™' = 4.0 GeV , V"^^ = 0.40 , d = -8.78 , Md = 0.97 GeV . 

Note that in this fit the large value of M™' is correlated with the small value of r^'^^ re- 
minding us that the functional t-dependence is not very well constrained from present data. 
The results of our two "dispersion-relation" fits and three hybrid model fits are available 
as a computer program providing the four-fold cross section of polarized lepton scattering 
on unpolarized proton for a given kinematics, see http://calculon.phy.hr/gpd/, Unlike 
"dispersion-relation" fits, the hybrid model fits, where LO evolution of sea quark and gluon 
GPDs has been taken into account, are suitable for estimates in the small xbj region. 

In Fig. [5] we confront our fit results (x^/d.o.f. ~ 1 w.r.t. the employed data sets) to 
experimental data: KM09a (squares), KM09h (circles), and the hybrid model fit KMlOb 
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Figure 6: SmH/vr obtained from different strategies: our DVCS fits [dashed (solid) curve excludes 
(includes) Hall A data from "dispersion relation" KM09a {KM09b) [S3] and hybrid KMlOb (dash- 
dotted) models], GK 07 model from DVEM (dotted), seven-fold CFF fit [92 [93] with boundary 
conditions (squares), TL, "H CFF fit [90 (diamonds), smeared conformal partial wave model fit [94] 
within H GPD (circles). The triangles result from our neural network fit, cf. Fig. [7| (left). 

(triangles-up) in which we now utilized the improved formulae set [87| and the Kelly form 
factor parameterization [ID] . Here we also include the predictions from the GK07 model 
(triangles-down), where we adopt the hypothesis of H dominance. Qualitatively, these 
predictions are consistent with a vggEI code estimate, which tends to over-estimate the 
BSAs isni [86] and describes the EGAs from HERMES rather well without L>-term [79] , 
This is perhaps not astonishing, since the employed H GPD model relies on RDDA, too. 
We would hke to emphasize that at LO the GK07 model is in reasonable agreement with 
the HI and ZEUS DVGS data (x^/d.o.f. ~ 2), essentially thanks to the rather small and 
stable skewness ratio r^^^ of sea quarks, see dashed curve in the left panel of Fig. [3l 

Longitudinally polarized target data from GLAS [88] and HERMES [89] provide a handle 
on H [11], where mean values of GFF fits [90] in JLAB kinematics give two to three times 
bigger H contribution compared to our expectations (r^ ~ 1,6^ ~ 2). These findings are 
one to two standard deviations away from our big H ad hoc scenario of the KM09b fit, 
which is indeed disfavored by longitudinally polarized proton data. We like to add that 
with our present hybrid model a reasonable global fit, such as KM 10 above, is possible. In 
such a fit the Hall A data require a rather large pion pole contribution, inducing a large 
DVGS cross section contribution. Still we have not included the transversal target data 
from the HERMES collaboration [79j or the neutron data from Hall A [91j. 

So far we did not study model uncertainties or experimental error propagation, since 
both tasks might be rather intricate. To illuminate this, we compare in Fig.[6]our results for 
9m'H(xB, i)/7r versus xb at i = —0.28 GeV^ (left) and for Hall A kinematics xb=0.36 versus 
—t (right) with results that do provide error estimates. The squares arise from constrained 
least squares fits \92\ [93] at given kinematic means of HERMES and JLAB measurements 
on unpolarized proton, where the imaginary and real parts of twist-two GFFs are taken as 
parameters. Note that Qmf and the other remaining eight GFFs are set to zero, however, all 
available observables, even those which are dominated by these GFFs, have been employed. 
This might increase 'statistics', however, yields also a growth of systematic uncertainties. 

^ VGG refers to a computer code, originally written by M. Vanderhaeghen, P. Guichon, and M. Guidal. 
To our best knowledge the code for DVCS, presently used by experimentalists, employs a model that adopts 
RDDA [83]. 
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Figure 7: Neural network extraction of 5Re'H(xBj, i)/"" i^oni BCA [86] and BSA ^80^ data. 



The huge size of the errors mainly shows to hmited accuracy with which H can be extracted 
from unpolarized proton data alone [TT]. A pure H GPD model fit [9l] (circles) to JLAB 
data provides much smaller errors, arising from error propagation and some estimated model 
uncertainties. All three of our curves are compatible with the findings [92 ^ [93] and the H 
GPD model analysis [94J of CLAS data. However, for Hall A kinematics the deviation of the 
two predictions that are based on the H dominance hypothesis, see dashed curve and circles 
in the right panel, are obvious and are explained by our underestimation of cross section 
normalization by about 50%. Moreover, the quality of fit |94j x^/d.o.f. ~ 1.7, might provide 
another indication that CLAS and Hall A data are not compatible when this hypothesis 
is assumed, see, e.g., the two rightmost circles in the left panel for CLAS (xbj = 0.34, 
t = -0.3GeV2, Q2 ^ 2.3GeV2) and Hah A (xbj = 0.36, t = -0.28GeV2, Q2 ^ 2.3GeV2). 
The pure Ti and Ti CFF fit [90] (diamonds), including longitudinally polarized target data, 
is within error bars inconsistent with the H dominated scenario [94] (circles), however, 
(accidentally) reproduces our dashed curve. 

Another source of uncertainties are twist-three contributions and perhaps also gluon 
transversity related contributions, which might be strongly affected by twist-four effects [95]. 
The Hall A cross section measurements allow us to have a closer look to the spectrum of 
harmonics and we emphasize that the second harmonics in Hall A data, i.e., effective twist- 
three contributions, are tiny or hard to separate from noise. Such contributions are smalU 
in HERMES kinematics, too, where the constant A^^, appearing at twist-three level, is a 

twist-two dominated quantity that, as expected [H], turns out to be correlated with A^^. 
However, we emphasize that even small twist-three effects, e.g., ~ 2%, might induce a larger 
uncertainty in a twist-two related quantity, e.g., ~ 10%. 

All this exemplifies that within (strong) assumptions and the present set of measure- 
ments the propagated experimental errors cannot be taken as an estimate of GPD uncer- 
tainties. An error estimation in model fits might be based on twist-two sector projection 
technique |llj . boundaries for the unconstrained model degrees of freedom, and error prop- 
agation in the twist-two sector. Alternatively, neural networks, already successfully used 
for PDF fits [96], may be an ideal tool to extract CFFs or GPDs. We present in Fig. [7] a 
first example in which, within //-dominance hypothesis, Ti is extracted using a procedure 
similar to the one of Ref . [97J . Here 50 feed- forward neural nets with two hidden layers were 
trained using HERMES BCA [86] and CLAS BSA [80j data. Hence, only the experimental 
errors were propagated, which in absence of a model hypothesis get large for the t — )> 

^Except for 3x2 beam spin asymmetry data points at largest —t, xb, and Q^, respectively [86| . 
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Figure 8: KMlOb model estimate for the DVCS beam spin asymmetry (PT|) with a proton (sohd) and 
neutron (dashed) target. Left panel: Abs versus (p for En = 250 GeV, = 5 GeV, xbj = 5 x 10"'^, 
= 10 GeV^, and t = —0.2 GeV^. Right panel: Amplitude Agg of the first harmonic versus xbj 
at t = -0.2 GeV^ for small xbj (thin) [E^ = 30 GeV, Ep = 360 GeV, ^ 4 GeV^] and large sej 
(thick) [Ee = 5 GeV, = 150 GeV, = 50 GeV^] kinematics. 



extrapolation. 

4 Potential of an electron-ion colider 

A high luminosity machine in the collider mode with polarized electron and proton or ion 
beams would be an ideal instrument to quantify QCD phenomena. It is expected that such 
a machine, combined with designated detectors, would allow for precise measurements of 
exclusive channels. Besides hard exclusive vector meson and photon electroproduction, one 
might address the behavior of parity-odd GPDs Ti, related to polarized PDFs, and £ via the 
exclusive production of pions even in the small x region. It is obvious from what was said 
above that an access of GPDs requires a large data set with small errors. In the following 
we would like to illustrate the potential of such a machine for DVCS studies, where we also 
address the GPD deconvolution problem. 

Let us remind that already the isolation of CFFs is rather intricate. For a spin-1/2 target 
we have four twist-two, four twist-three, and four gluon transversity related complex valued 
CFFs. The photon helicity non-flip amplitudes are dominated by twist-two CFFs, the 
transverse-longitudinal flip amplitudes by twist-three effects, and the transverse-transverse 
flip ones by gluon transversity. Hence, the first, the second, and the third harmonics 
w.r.t. the azimuthal angle of the interference term are twist-two, twist-three, and gluon 
transversity dominated. In an ideal experiment, assuming that transverse photon helicity 
flip effects are negligible, cross section measurements would allow to separate the sixteen 
quantities that are then given in terms of twist-two and twist-three CFFs. The reader might 
find a more detailed discussion, based on a 1/Q expansion, in Ref. [H]. We add that the 
definition of CFFs is convention-dependent. 

In a twist-two analyzes on unpolarized, longitudinally and transversally polarized pro- 
tons one might be able to disentangle the four different twist-two CFFs via the measurement 
of single beam and target spin asymmetries. In Fig. [8] we illustrate that for a proton target 
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Figure 9: DVCS longitudinal target spin asymmetry versus (j) for KM09a (dashed), KM09b (solid), 
and KMlOb hybrid (dash-dotted) models at ^ 5 GeV, t = -0.2 GeV^ ^ 4 GeV^ within 
Ep = 150 GeV, XBi = 0.1 (left) and Ep = 350 GeV, xb; = 0.01 (right). 



(solid) the beam spin asymmetry (j2ip 

A\^l^y\Fi{t)H{ttt,Q^) 



4M2 



F2{t)E{C,C,t, Q2) + 



(24) 



might be rather sizeable over a large kinematical region in which the lepton energy loss y 
is not too small. Here the helicity conserved CFF H is the dominant contribution, while 
£ appears with a kinematic suppression factor t/AMf^, induced by the helicity flip. For a 
neutron target the T-L contribution is suppressed by the accompanying Dirac form factor 
{F^{t = 0) = 0) and so one becomes sensitive to the CFF £. Unfortunately, one also has 
to worry about other non-dominant CFF contributions, indicated by the ellipsis. Note that 
the asymmetry for neutron (dashed) might be underestimated, since we set in our model 
E{x,x,t, Q^) to zero. 

For a longitudinally polarized target the asymmetry 
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(25) 



is sensitive to the GPD H, where and other GPDs might contribute to some extent. 
Naively, one would expect that this asymmetry vanishes in the small xb] region and might 
be sizeable at xbj ~ 0.1, see left panel of Fig. [9l Not much is known about the small x 
behavior of H and it might be even accessible at smaller values of xbj , as illustrated by the 
KM09h model with its big H contribution (solid, right panel). For a neutron target the 
asymmetry becomes sensitive to the ^E GPD. Note that here the factor ^ is annulled by a 
conventional 1/^ factor in the definition of E GPD. 

Finally, we emphasize that a single spin asymmetry measurement with a transversally 
polarized target provides another handle on the helicity- flip GPDs E and E GPDs. If the 
target spin is perpendicular to the reaction plane, the asymmetry 

^ [F^mi.^, e, Q') - F^{t)E{i, ^, t, q2) + . . . ] ^ (26) 
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is dominated by a linear combination of the GPDs H and E. In the case that the target 
spin is aligned with the reaction plane the asymmetry 

t 



A 



m 

■TS 



OC 



4M2 



F2{t)H{i, t, Of) - Fi{t)^E{C, e, t, Q') + 



(27) 



15 



0.4 



0.2 



xsj=5 10 
-1=0.2 GeV^ 
Q^=10 GeV^ 



J? 0.0 



-0.2 




-0.4 



Figure 10: KMlOh model estimate for tlie DVCS beam clrarge asymmetry witli a proton (solid) 
and neutron (dashed) target. Left panel: ^bc versus (j) for Em = 250 GeV, = 5 GeV, x-gj — 
5 X 10"^ Q2 ^ iQ QeV^, and t = -0.2 GeV^ Right panel: Amplitude A^^^l of the first harmonic 
versus xb; at t = -0.2 GeV^ for smaU xbj (thin) [E^ = 30 GeV, Ep = 360 GeV, = A GeV^] and 
large sbj (thick) [E^ = 5 GeV, Ep = 150 GeV, = 50 QeV^] kinematics. 



is dominated by a linear combination of the GPDs H and E. Unfortunately, these asymme- 
tries are kinematically suppressed by the factor t / 4M^ and for a neutron target in addition 
by the Dirac form factor Fi(t). 

Although the given formulae (|21H27p are rather crude, they illustrate that a measure- 
ment of single spin asymmetries would allow to access the imaginary part of the four twist- 
two related CFFs, however, the normalization of these asymmetries depends to some extent 
also on the real part of the twist-two related CFFs and the remaining eight ones. Measure- 
ments of cross section differences would allow to eliminate the normalization uncertainty, 
and combined with a harmonic analysis one can separate to some extent twist-two, twist- 
three, and gluon transversity contributions. However, also then the extracted harmonics 
might be contaminated by DVCS cross section contributions, which are bilinear in the 
CFFs. To get rid of these admixtures, one needs cross section measurements with a positron 
beam. Forming differences and sums of cross section measurements with both kinds of lep- 
tons, allows to extract the pure interference and DVCS squared terms and might allow so 
to quantify twist-three effects. Existing data indicate that these are small, as is expected 
based on kinematic factors. However, even obtaining only an upper limit is important for 
determination of the systematic uncertainties of twist-two CFFs. 

We also emphasize that having both kinds of lepton beams available allows to measure 
the real part of CFFs. In Fig. [10] we show the beam charge asymmetry ()22p . 



(28) 



for an unpolarized target, which is expected to be sizeable. For a proton target this asym- 
metry should possess in the transition from the valence to the sea region a node (thick solid 
curve, right panel). In our parameterization the real part of the £ CFF is determined by 
the P subtraction term, which induces even for neutron target a sizeable asymmetry (thick 
dashed curve, right panel). 

The large kinematical coverage of the proposed high-luminosity EIC raises the question: 
Can one utilize evolution, even at moderate xbj values, to access GPDs away from their 
cross-over line? Similarly as it has been done for the small xbj region, we use the Mellin- 
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Figure 11: The upper left panel shows the valence- like contribution (|T9)) to the CFF % extracted 
with a "dispersion-relation" fit KMOQa from fixed target measurements (dotted) aX t — —0.2 GeV^ 
and — 2 GeV^ versus xb] together with various models. The lower left panel shows the cor- 
responding GPD models xH{x,ri,t, Q^) together with a minimalist GPD parameterization (dotted 
curve) versus a; at 77 = 0.2, t — —0.2 GeV^, and — 2 GeV^. The corresponding quantities at 
Q"^ — 50 GeV^ are displayed in the right panels. 



Barnes integral technique to address the problem. Taking different non-leading S0(3) par- 
tial waves in the ansatz for the conformal moments (|14|15p . we build three different GPD 
models for valence quarks that provide almost identical CFFs, see upper left panel in Fig.llll 
They are compatible with both, PDF and form factor parameterizations, see Fig. [H and 
the outcome ()19p from the "dispersion-relation" fit KM09a (dotted curve). We note that 
the different model behavior at large xbj results only in a smaller discrepancy for the real 
part of the CFF in the kinematics of interest. In the lower left panel of Fig. [11] we illustrate 
that for fixed 77 the x-shape of the three GPD models looks quite differently. Compared to 
the minimalist model (dotted curve), a model with a negative next-to- leading partial wave 
(solid) decreases the GPD size on the cross-over line r] = x and generates an oscillating be- 
havior in the central region. Model with an alternating-sign S0(3) partial wave expansion 
(dash-dotted) possesses more pronounced oscillation effects in the central region or even 
nodes. In the third model (dashed curve), the reduction on the cross-over line is reached 
within a next-to-next leading S0(3) partial wave. Note that the GPDs in the region r] x 
are governed by the x-behavior of the PDF analogues. In the right panels we demonstrate 
that for a large lever arm in Q'^ evolution effects e.g., at = 50 GeV^ are important in 
the valence quark region. However, for CFFs (upper right panel) the discriminating power 
of evolution effects remains moderate even if the GPD shapes look rather different. 
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5 Conclusions and summary 



With all the theoretical tools sketched above plus those which are presently under develop- 
ment it is clear that our understanding of hadron structure will be revolutionized once most 
of the diverse asymmetries can be measured with percent or permille precision (depending 
on the observable). A high-luminosity EIC, as it is proposed at Brookhaven National Lab., 
would allow to do so. 

Let us summarize. At present, the first steps have been undertaken to access GPDs from 
experimental data in the small a:;Bj region with models and in the fixed target kinematics 
within different strategies, providing some insight into the GPD H. In particular for hard 
exclusive photon electroproduction in the fixed target kinematics, model fits to leading order 
accuracy in a^, relying on the scaling hypothesis, are compatible with least square CFF fits 
and with a first result from neural networks, assuming H dominance. The large uncertainties 
here in extracting CFFs, including "H, are mainly related to the lack of experimental data. 
Thus, not only the extraction of the very desired £, playing an important role in the "spin- 
puzzle", but also of other CFFs, requires a comprehensive measurement of all possible 
observablcs in dedicated experiments. The comparison of H GPDs accessed from DVCS at 
leading order with model fits and from hard exclusive meson production in the hand-bag 
approach within Radyushkin's double distribution ansatz shows that in the valence region 
the extracted quark GPDs are somewhat different while at small x they are becoming 
compatible. The main difference lies in the gluonic sector and is induced by utilizing PDF 
parameterizations that where extracted in leading and ncxt-to-lcading order from inclusive 
measurements. A more appropriate analysis of these processes requires the inclusion of 
radiative corrections in a global fitting procedure, which is in progress. We should also 
mention here that hard exclusive processes with nuclei, which at present are not extensively 
studied, opens a new window for a partonic view on their content. 

Imaging the partonic content of the nucleon and the phenomenological access to the 
proton spin sum rule from hard exclusive processes can only be reached through proper 
understanding of GPD models. It might be pointed out here that GPDs can also be for- 
mulated in terms of an effective nucleon (light-cone) wave functions, which links them to 
transverse momentum dependent parton distributions. A whole framework, consisting of 
perturbative QCD, lattice simulations, and dynamical modeling, is available to reveal GPDs 
and to access the nucleon wave function. Such a unifying description might be considered 
as the primary goal, which quantifies the partonic picture. Although such a task looks 
rather straightforward, much effort is needed on the theoretical, phenomenological, and 
experimental side, where experimental data with small uncertainties play the key role. A 
high-luminosity EIC is an ideal machine that would cover a wide kinematical range, would 
complement the planned fixed target experiments at JLAB@11 GcV, and has, besides new 
measurements, great potential to significantly improve existing data sets. 
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